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1. Introduction 

Using powerful simulation software has the advantage of including many effects and potentially 
results in high accuracy. On the other hand it does not always give the insight an equation can 
deliver. 

This paper originates from the necessity to understand the Pulse Height Spectra (PHS) given by 
solid neutron converters employed in thermal neutron detectors as in [|]], [Q], and from the 

investigation over such a detectors' efficiency optimization. Many efforts have been recently made 
in order to address the ^He shortage problem. The development of new technologies in neutron 
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detection is important for both national security [||] and for scientific research Examples of 
application can be found in [0], and [g]. 

When a neutron is converted in a gaseous medium, such as a ^He detector, the neutron capture 
reaction fragments ionize the gas directly and the only energy loss is due to the wall effect. As 
a result, such detectors show a very good gamma-rays to neutron discrimination because gamma- 
rays release only a small part of their energy in the gas volume and consequently neutron events 
and gamma-rays events are easily distinguishable on the PHS. 

On the other hand, when dealing with hybrid detectors as in [p, where the neutron converter is solid 
and the detection region is gaseous, the gamma-ray to neutron discrimination can be an issue [p^], 
[^]. Indeed once a neutron is absorbed in the solid converter, it gives rise to charged fragments 
which have to travel across part of the converter layer itself before reaching the gas volume to 
originate a detectable signal. As a result, those fragments can release only a part of their energy 
in the gas volume. The neutron PHS can thus have important low energy contributions, therefore 
gamma-ray and neutron events are not well separated just in energy. 

In this paper we want to give a comprehension of the important aspects of the PHS by adopting a 
simple theoretical model for solid neutron converters. We will show good agreement of the model 
with the measurements obtained on a ^''B-based detector. 

In the same way the analytical model can help us optimize the efficiency for single and multi-layer 
detectors in different circumstances of incidence angle and neutron wavelength distribution. 



The model we use is the same as implicitly used in many papers such as [ ]11[ ] or [|12|]. It makes the 
following simplifying assumptions: 

• the tracks of the emitted particles are straight lines emitted back-to-back and distributed 
isotropically; 

• the energy loss is deterministic and given by the Bragg curves without fluctuations; 

• the energy deposited is proportional to the charge collected without fluctuations. 

Referring to Figure we talk about a back-scattering layer when neutrons are incident from the 
gas-converter interface and the escaping particles are emitted backwards into the gas volume; we 
call it a transmission layer when neutrons are incident from the substrate-converter interface and 
the escaping fragments are emitted in the forward direction in the sensitive volume. We consider 
a neutron to be converted at certain depth (x for back-scattering or d — y for transmission) in the 
converter layer and its conversion yields two charged particles emitted back-to-back. 

2. Double layer 

We put a double coated blade in a gas detection volume. A blade consists of a substrate holding 
two converter layers, one in back-scattering mode and one in transmission mode. 



Starting from the analytical formulae derived in Ql IQ we are going to derive properties that can 
help to optimize the efficiency in the case of a monochromatic neutron beam and in the case of a 
distribution of neutron wavelengths. 

By denoting with dss the thickness of the coating for the back-scattering layer and with dj the 
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Figure 1. Variables definition for a back-scattering and transmission layer calculations. 



transmission layer thickness, the efficiency of the whole blade is: 



e{dBs,dT) = £bs (dBs) + e ^'^"^ ■ £t [dr) 



(2.1) 



where £^5 {dss) and ej {dr) are the efficiency for a single coating calculated as shown in appendix 
^ from [jn]], and £ is defined in Equation |2.3| . The relation Ve{dBs,dT) = determines the two 
optimal layer thicknesses. 

In order to keep calculations simple, we consider only two neutron capture fragments yielded by 
the reaction. This approximation will not affect the meaning of the conclusion. In the case of ^Li 
Equation 01 is exact, for ^^B the expression O] should ideally be replaced by e{dBs,dT) = 0.94- 



+ f'-^-^si- • £0-94 [dT)) + 0.06 • (£*]j'^ {dBs) + e-^'^^' ■ £?"^ {dT)); where £*'-94 means the 
efficiency calculated for the 94% branching ratio reaction with the right effective particle ranges. 
We will limit us to the 94% contribution as if it were 100%. R\ and R2, with {R2 < Ri), are the 
two ranges of the two neutron capture fragments. In case of ^"B4C the two 94% branching ratio 
reaction particle ranges are. Ri =3 pLm (a-particle) and /?2 = 1.3 pLm CU), when a 100 KeV energy 



threshold is applied (as defined the minimum detectable energy in []11[]). 

As £bs{'^) and £7- (<i) have different analytical expressions according to whether d < R2 < R\, 
R2 < d < Ri or R2 < Ri < d we need to consider 9 regions to calculate Ve{dBs,dT) as shown in 
Figure ^ If we were to include the four different reaction fragments we would have to consider 
25 domain partitions. We will see later that the important regions, concerning the optimization 
process, are the regions square 11 (where dBs < ^2 < ^1 and dj < R2 < Ri) and square 22 (where 
R2<dT< Ru R2 < dBs < Ri)- 

In order to consider a non-orthogonal incidence of neutrons on the layers, it is sufficient to replace 
£ with -^^0^, where 6 is the angle between the neutron beam and the layer surface (see Figure [I]). 
This is valid for both the single blade case and for a multi-layer detector. The demonstration can 



be found in [ 11]. 
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Figure 2. Domain for the blade efficiency function. The domain is divided into 9 partitions according to the neutron 
capture fragment ranges. In red and blue are the domain partitions where the efficiency directional derivative along 
the unity vector u is respectively always positive or always negative. The red dotted line represents the case when the 
substrate effect is not negligible and the maximum efficiency can not be attained on the domain bisector. 



2.1 Monochromatic double layer optimization 

In the domain region called square 11 the efficiency turns out to be: 

£11 idBS,dT) =A- (l-e-^'^'^A + des ■ C ■ e-^'^'' + e'"''^'^' (b- (l -e"^^^ 



C-di 



(2.2) 



Where we have called: 

1 ' ' 



B 



ILRi ILR2 



C 



1 1 



E = ?i-a(A) 
(2.3) 



where n is the number density of the converter layer and CJ(A) its neutron absorption cross-section. 
By calculating Vei 1(^55,^7-) = we obtain the result that dgs = dj and 



(2.4) 



We repeat the procedure for the square 22 in which the efficiency is: 

£22{dBS,dT) =A + §zr7~^T^" 

+e-^'^ ■ {-B ■ e-^^ + i (1 + i - t ) + 
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We obtain again dgs = dj and 



T-BS 



1 1 ( ^ 



(2.6) 



Naturally each result of Equations 2.4 and is useful only if it gives a value that falls inside the 
region it has been calculated for. 



The points defined by Equations ^ and define a maximum of the efficiency function in the 
regions square 11 or square 22 because the Hessian matrix in those points has a positive determi- 
nant and -^-S- is negative. It is easy to demonstrate there are no extreme points outside the domain 
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Figure 3. Efficiency plot for a double coated substrate with '"i?4C at = 90° incidence at 1.8A (left) and 20A (right). 
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Figure 4. Efficiency plot for a double coated substrate with '"i?4C at = 10° incidence at 1.8A (left) and 20A (right). 

regions where either dss > Ri or dj > Ri, i.e. square jk with j = 3 or A; = 3 or both. This outcome 
is also intuitive. In back-scattering mode when the converter thickness becomes thicker than the 
longest particle range there is no gain in efficiency by adding more converter material. In the 
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transmission case, increasing the thickness above Ri will add material that can only absorb neu- 
trons without any particle escaping. 

For the cases of squares 12 and 21, we obtain that Ve{dBs,dT) = has no solution; thus the effi- 
ciency maximum can never fall in these domain regions. 

In Figures |3| and ^ the efficiency for four different cases is plotted. The red circle identifies the 
point of maximum efficiency calculated by using Equations or it stands out immediately 
that even though the efficiency function is not symmetric relatively to the domain bisector (drawn 
in black) the maximum nevertheless always lies on it. 

This is a very important result because the sputtering deposition method /|7j|/ coats both sides of 
the substrate with the same thickness of converter material and it is also suited to make optimized 
blades. 

2.2 Effect of the substrate 



If we consider the neutron loss due to the substrate, the Equation |2JJ has to be modified as follows: 



e{dBs,dT) = Sbs {dss) + 5(A) • e'^'^'^' ■ Et {dr) (2.7) 

Where 5(A) = and 'Lgut and d^ub are the macroscopic cross-section and the thickness 

of the substrate. 



If we optimize the layer thicknesses, we find the same result of Equations and for the 
transmission layer thickness dj but, on the other hand, the back-scattering layer thickness does not 
equal the transmission layer thickness anymore. It becomes: 

^ ■ {\ — d) + d -dr for square 11 

Ri-{l-25) + 5-(l + ^) -dr for square 21 (2.8) 



Ri • {I — 5) + 5 ■ dr for square 22 

The maximum efficiency can now also lie in square 21 but not in square 12, because 5 G [0, 1]. 
The dotted line in squares 11, 21 and 22, in Figure ^ are Equations |2.8[ The slope of the line in 
squares 11 and 22 is equal to 5(A) < 1. The maximum efficiency, when 5(A) is not negligible, 
now lies on the dotted line just above the optimum without substrate effect. 

From Equation we can observe that when 5 is close to zero, i.e. the substrate is very opaque to 
neutrons, the thickness of the back-scattering layer tends to the value R\. On the other hand, when 
5 is close to one, dss tends to dj- The factor 5(A) is usually very close to one for many materials 
that serve as substrate. We define the relative variation between dss and dj as = ^^^^'^^^ ; in 

Table |l| they are listed for a 0.5mm and 3 mm Aluminium substrate of density p = 1.1 gjcnr' at 
1.8A. In case the substrate is inclined under an angle of 10° a substrate of 0.5 mm looks like a 
substrate of about 3 mm. We consider a neutron to be lost when it is either scattered or absorbed, 
therefore, the cross-section used [|l4|] is: Om = offiX) + ajf " = 0.2Z?(at 1.8A)+1.5Z7 = \nb. 



2.3 Double layer for a distribution of neutron wavelengths 

The result of having the same optimal coating thickness for each side of a blade is demonstrated 
for monochromatic neutrons and we want to prove it now for a more general case when the neutron 
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Table 1. Neutron loss factor S for an Aluminium substrate and relative difference between the two coating thicknesses 
held by the substrate for 1.8 A. 



dsuh{mm) 


5 (1.8A) 


\d (1.8A) 


0.5 

3 


0.995 
0.970 


0.0004 
0.0026 



beam is a distribution of wavelengths and when the substrate effect can be neglected (5(A) « 1). 
We will prove a property that will turn out to be useful. We will show that the directional derivative 
of E{dBs,dT) along the unit vector u = is positive until the bisector dss = dr and it 

changes sign only there. This vector identifies the orthogonal direction to the bisector (see Figure 

§. 

In the square 11: 

DuE = ux We{dBs,dT) = -j= ■ e'^'^' (dgs - dj) (2.9) 

v2 



In the square 22: 



Dae = ■ e-^'^'^ (des - dr) (2.10) 



Which are both positive above the bisector and negative below. In the other domain regions the 
demonstration is equivalent. E.g. in the square 12: 

2^2 V ^2 Ri RiJ 



Which is strictly negative in the square 12 where dss < R2 < Ri and R2 < dj < Ry except in the 
comer, on the bisector, where dss = dT =R2- 
The following theorem is therefore proved. 

Theorem 2.1 The efficiency function defined by the Equation ^ is strictly monotone in the two 
half-domains identified by the bisector dss = dj (see Figure 

In a general instrument design one can be interested in having a detector response for a whole 
range of A. E.g. an elastic instrument can work a certain time at one wavelength and another 
time at another wavelength. In a Time-Of-Flight instrument one can be interested in having a 
sensitivity to neutrons of a certain energy range including or excluding the elastic peak. One can 
define a normahzed weight function w (A) (/q^°°w (A) dX = \) that signifies how much that neutron 
wavelength is important compared to others. I.e. the price we want to spend in a neutron scattering 
instrument to be able to detect a neutron energy with respect to an other one. We can also consider 
an incident beam of neutrons, whose wavelength distribution is w(A), and we want to maximize 
the efficiency given this distribution. 

The efficiency for a blade exposed to a neutron flux which shows this distribution is: 

/.+00 

ew{dBs,dT)= w{X)e{dBs,dT,X)dX (2.12) 
Jo 
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where e{dBs,dT,X) is the efficiency in Equation 2A 



In order to optimize this efficiency its gradient relative to dss and dj has to be calculated: 

Ve^{dBs,dT)= w{X)Ve{dBs,dT,X)dX (2.13) 

^0 

Both gradient components have to cancel out {-j^ = ^ = 0), this leads to D^Ew = 0. E.g. in 
square 11: J^°°w{X) (^e^^^'^^ {dj —dBs)CL) dX =0. As a result, in order for the efficiency to 
attain a maximum, it is necessary (but not sufficient) that its directional derivative along the unity 
vector M = (1, — 1): 

Dae^ = / w{X)Dae{dBs,dT,X)dX (2.14) 

JO 

equals zero. 

For a general family of functions f{dBs,dT,X) for which the maximum always lies on the domain 
bisector it is not true that the function defined by their positively weighted linear combination 
must have the maximum on dBs = dj, because in general Vf{dBs,dT,X) can be positive, null or 
negative, thus there are many ways to accomplish Vfw{dBSidT) = 0. However, thanks to Theorem 
|2.l| , DiiE = is satisfied only on the bisector. Below the bisector, this is always negative, as it is a 
positively weighted integral of negative values; similarly, above the bisector, this is always positive. 
Hence, the maximum can only be attained on the bisector. 

The gradient can hence be replaced by and the function maximum has to be searched on the 
bisector, therefore: 

VeUdT)= r w{X)^e{dBs = dT,dT,X)dX (2.15) 
JO odj 

In the end, the same layer thickness for both sides of a blade has to be chosen in order to maximize 
its efficiency, even if it is exposed to neutrons belonging to a general wavelength distribution w{X). 
The integration over X can be alternatively executed in the variable £; indeed £ is just a linear func- 
tion in X because Oahs is proportional to X in the thermal neutron region. Moreover, as indicated 
previously, £ is also a function of d and this is the only appearance of 6 in the efficiency function. 
Hence, we can just as well consider a weighting in X and Q which results in just a weighting func- 
tion over £. In other words, all the results that have been derived for a wavelength distribution also 
hold for an angular distribution or both. 

2.3.1 Flat neutron wavelength distribution example 

As a simple example we take a flat distribution between two wavelengths X\ and X2 defined as 
follows: 

^(^) = ^— ^ (2.16) 
A2 — Ai 



In the square 11 we obtain: 
_d 



^ E{dBs = dT,dT,X)= ^e-^^ (^Ble^^^ - c) (2. 17) 
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We call El = and £2 = ^(■^2)- We recall that A and B are function of 



1 



, . ( 2Ce+^^-C-r--^- ) 



1:2 



(2.18) 



In the same way the solution in the square 22 can be determined. 



de 
ddr 



e-^^ i e-^^ ( 2BL- 



+1.R2 
R2 ) R\ 



(2.19) 



By integrating we finally obtain: 
1 



V£w{dT 



(S2 



C-Z 



1 







(2.20) 



dj V ^1 2(3fr /?2(/?2-2(ir). 

The solution of Equations [2.18 and 2.20 gives the optimum value for the thickness of the two 
converter layers in the region of the domain called square 11 and square 22 respectively for a 
uniform neutron wavelength distribution between Ai and A2. E.g. for a uniform neutron wavelength 
distribution between lA and 20A the optimal thickness of coatings on both sides of the blade is 
1 nm. 



3. The multi-layer detector design 



12 3 4 5 6 
1 2 3 



substrate 



coatings 



Figure 5. A multi-layer detector schematic. N/, blades (holding N = 2-Ni, converter layers) are placed in cascade, 
alternating with detection regions. 

In a detector Uke that presented in [jl]], [|l5]], or in [16], all the substrates have the same coating 
thickness. One can ask if it is possible to optimize the coating thicknesses for each layer in order to 
gain in efficiency. This is also applicable to neutron detectors which use solid converters coupled 



with GEMs []17[]. In a multi-layer detector (see Figure g), composed by N layers ov Nh= 2 blades, 
the efficiency can be written as follows: 



etot{Nb) = ei{dBsi,dTi) + £ eiidBsk,dTk) -e (^.'AV'(^«s,+dr,)) 

k=2 



■Z 



(3.1) 
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Where Si{dBsk,dTk) = ^ik represents the efficiency for a single blade already defined by the Equa- 
tion dssk and djk are the coating thicknesses of the k — th blade. 

If the detector is assembled with blades of the same thickness, i.e. dssk = djk = d, \/k = l,2,... ,Nb, 



Equation 3.1 can be simplified as follows: 



etot{Nh) = Slid) -J^e 



-2{k-\)d-Z 



ey{d)- 



-IdlNb 



k=\ 



1 



-2dl. 



(3.2) 



Therefore, -ff- = has to be calculated in order to optimize such a detector containing blades of 
same coating thicknesses. 



In the case of a distribution of wavelengths, defined by w(A), Equation ^Jjhas to be integrated 
over X as already shown in Section 2.3 . 

eZANh,dBs,d~T) = w{X)etot{Nb,X)dX (3.3) 
Jo 

The efficiency will be function of N = 2 ■ Nb variables; which can be denoted using the compact 
vectorial notation by the two vectors Jbs and dj of Nb components each. 

Both for monochromatic mode and for a distribution the detector optimization implies the calcula- 
tion of the A/^-dimensional gradient VStot, because in the case of a distribution the gradient can be 
carried inside the integral over X. We will use the shorthand Sik = £i {dBskidrk)- 
The k — th A/^-dimensional gradient component for back-scattering mode results to be: 



d£t, 
dd, 



BSk 



dssk 



ddBSI ^ 



ddssk 



l!^T\dBSj+dr,))-Z 



i:,'=i''(^(M,+rfr,))-2: 



if A:= 1 

{P^'^ids.^dr.))-Z if 1<^<A^ 

if k = Nb 

(3.4) 



Equivalently we find the same formal expression for the k — th component of the gradient with 



respect to the transmission variable {dj)', we can substitute ddBsk with ddrk in Equation 3A which 
are for the rest entirely the same. This implies that if we put identical conditions on 



dds: 



and on 



de, 



this will result in 

odjk adss 



deik — deik 



ddj 



. Equivalently as already found in Section 2.3, from Equation 3.4 



and the one for the transmission variable we finally obtain {\/k = 1,2, .. . jNb)'- 



deik 



ddssk ouTk 
deik — deik 
ddgsk ddrk 



if k<Nb 
if k = Nb 



(3.5) 



Both for the monochromatic case and in the case of a distribution of wavelengths, the condition in 



Equation ^ has to be satisfied. Using Theorem |2.1| , the maximum efficiency can only be found, 
again, on the bisector. 

In a multi-layer detector, which has to be optimized for any distribution of neutron wavelengths or 
for a single wavelength, all the blades have to hold two layers of the same thickness. Naturally, 
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thicknesses of different blades can be distinct. 

Thanks to this property, we can denote with dk the common thickness of the two layers held by the 
k — th blade {dssk = djk = dk), furthermore, Equation 3. 1 can be simplified as follows: 



Etot {N, d) = £1 ( ) + £ £i {dk) ■ e 
k=2 



-2 Ei^i 



(3.6) 



where d is the vector of components dkfovk= 1,2,. .. ,Nb. 

Optimizing a detector for a single neutron wavelength or for a distribution is different; the equation 
Vetot = in one case and Ve/^, = in the other represent a A//, -dimensional system of equations 
in Nb unknowns because of the simplification of having the same back-scattering and transmission 
layer thickness on the blades. By expanding the Equation |3^ we obtain: 



etotiNb) = ei{di)+e 



Slid: 



e\{d^) + ... 



+ e 



-2dil. 



(3.7) 



ei(Ji)+e-2^'^- ei{d2)+e 



^-2dil. 



We notice that the variable Ja?,^ appears only once, this means that, in the case of a single wave- 
length, its value can be determined without taking the others into account. Continuing the reasoning 
we see that the system of equations is upper triangular. In the monochromatic case we can opti- 
mize the detector starting from the last blade and going backward till the first. This is not true for 



the distribution case in which the gradient of Equation 3.7 is in addition integrated over A, thus 
all the blades have to be taken into account simultaneously in the optimization process. In the 
monochromatic case, we can start by fixing the last blade coating thickness because any change on 
the previous will only affect the number of neutrons that reach the last blade, and we require the 
last blade to be as efficient as possible for that kind of neutron. As the layer thickness optimum of 
each blade does not depend on the previous ones, the system is triangular. On the other hand, in the 
case of a wavelength distribution, any change on the previous blades will change the actual distri- 
bution of wavelengths the last blade experiences. Thus, the neutron distribution a blade has to be 
optimized for depends on all the previous blade coatings. In this case, the system is not triangular. 

3.1 Monochromatic multi-layer detector optimization 

In order to optimize the layers in multi-layer detectors for a given neutron wavelength, we can 
maximize the last layer efficiency and, then, go backward until the first layer. Formally, from 



Equation |3.7| , we obtain an iterative structure: 



fk 



' £i {dk) + e-^''^^ -ak+i ifk< Nb 



ei{dk 



if k = Nh 



(3.8) 



Uk+i is a fixed number, independent from dk, and represents the cumulative efficiency of the detec- 
tor from the blade {k+l) — th to the end. 



dh^filrMdk) 
ddk \^e,{dk) 



iik<Nb 
if k = Nb 



(3.9) 
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2g-£i (dk) are the derivatives in Equations 2.17 and 2.19| according to the domain partitions. In the 



domain region called square 11 as defined in Section ^, we obtain: 

dfk_i 2e-^^ {{B-au+i) - C) if ^ < A^z, 
ddk \2e-^''{B'Le-^i'-C) if k = Nh 



And in the square 22: 




if k = Nb 



(3.10) 



(3.11) 



Equations |3.10| and [3.11| have solutions similar to 2.4 and [2.6[ In the square 11 the solution is: 

c 



i . In , 



In the square 22: 



dT 



if k = Nb 



}..]„ El 1 

Ui I 2i?2i;(B-at+i)-<;+^2 

1 In ( ^ / 1 . ' 



if k<Nb 

if A: = A^/, 



(3.12) 



(3.13) 



The optimization method is a recursive procedure that employs the Equations p.l2| and |3.13| ; we 
start from the last blade, and we find its optimal thickness d'^'^' , afterwards we calculate Wa^,^ as the 
last layer efficiency using the optimal thickness found. Now we can calculate d'^'^J_^ from Equations 



3.12 or p. 13 and a^i^-i and so on until the first layer. 



CCk+l 



idr)-e 



('■-1 



-2 1- 



if k+l<Nb 
if k+l=Nh 



(3.14) 



3.1.1 Example of application 

We analyze a detector composed of 30 successive converter layers (15 blades) crossed by the neu- 
tron beam at 90° (like in Figure ^. We consider ^^840 (p = 2.24g/cm-') as converter; we neglect 
again the 6% branching ratio of neutron capture reaction. A 100 KeV energy threshold is ap- 
pUed and particle ranges turn out to he R\ =3 pLm (a-particle) and /?2 = 1-3 pim Chi), for the 94% 
branching ratio. 

Figures ^ and ^ show the optimization result for this multi-layer detector; for a monochromatic 



neutron beam of 1.8A and lOA. On the left, the optimal thickness given by either Equations 3.12 



or 3.13 is plotted in red for each blade; for comparison we use two similar detectors suitable for 



short and for long wavelengths in which the blades are holding l.ljxm and 0.5 jJLm thickness coat- 



ing. Those values have been obtained by optimizing the Equation the efficiency for a detector 
holding 15 blades of all equal thicknesses for 1.8A and for for lOA. The detector with 1.2/im 
coatings is very close to the one presented in On the right, in Figures || and ^ the efficiency 
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Figure 6. Thicknesses of the blade coatings (left) and their efficiency contribution (right), for a detector made up of 15 
blades of 1 .2 flm, 0.5 flm and for a detector optimized for 1 .8A. 
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Figure 7. Thicknesses of the blades coatings (left) and their efficiency contribution (right), for a detector made up of 
15 identical coating thickness blades of 1 .2 fim, 0.5 lim and for a detector optimized for lOA. 



Table 2. Efficiency for an optimized multi-layer detector and for a detector which contains 15 identical blades of 1 .2 /xm 
and 0.5 flm. 



wavelength (A) 


opt. detect. 


0.5 /im detect. 


1.2/im detect. 


1.8 


0.525 


0.388 


0.510 


10 


0.858 


0.831 


0.671 



contribution of each blade is plotted, again for an optimized detector for 1 .8 A and for an optimiza- 
tion done for lOA. The expression of the efficiency as a function of the detector depth is given by 
Equation 3.6 for each blade by fixing the index k. The whole detector efficiency is given in the 
end by summing all the blades' contributions. The whole detector efficiency is displayed in Table 
^for the detector of Figures ||and ^ By optimizing the detector for a given neutron wavelength we 
gain only about 2% efficiency which is equivalent to add more layers to the detectors optimized to 
hold identical blades. 
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In Figure ^ is shown the efficiency resulting from the monochromatic optimization process of the 
individual blade coatings and the optimization for a detector containing all identical blades (which 
thicknesses are shown on the right for each neutron wavelength). Neutrons hit the layers at 90° 
and five cases have been taken into account with an increasing number of layers. We notice that 
about for all neutron wavelengths the gain in optimizing the detector with different blades, let us to 
gain few percent in efficiency. The values in Table ^ are the values on the pink solid curve and the 
dashed one at 1.8 A and at lOA in Figure ||. 




Figure 8. Efficiency (left) and optimal thickness of the identical blades (right) as a function of neutron wavelength for 
a 2, 10, 20, 30 and 40 layers multi-layer detector. Solid lines indicate the optimized efficiency, for each wavelength, for 
a detector made up of blades of identical thicknesses; the dashed ones indicate the monochromatic optimization using 
different thicknesses. 



Still referring to Figure ||, we notice that a detector with 15 individually optimized blades (30 lay- 
ers) has about the same efficiency (above lOA) than a detector optimized to contain 20 blades (40 
layers) of equal thickness. On the other hand for short wavelengths the difference is not very sig- 
nificant. Moreover, there is also a trade off between the constraints of the detector construction and 
the complexity of the blade production. 



3.2 Multi-layer detector optimization for a distribution of neutron wavelengths 

In this case it is not possible to start the optimization from the last blade because the thicknesses 
of the previous layers will affect the neutron wavelength distribution reaching the deeper laying 
blades. We have in this case to optimize an A'^/, -dimensional function at once. Therefore, the Ny- 
dimensional equation Ve/^j = has to be solved: 

V,C=/ w{X)^dX=0 (3.15) 



is an expression similar to Equation 3.4 provided that we impose dssk = djk^k = 1,2, .. . ,Nb- 



In order to optimize a detector for a given neutron wavelength distribution w(A), the following 
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system of Nh equations in Nh unknown {dk) has to be solved: 



/o^"w(A) 



=0 
t/A =0 



t/A = if = 1 



(3.16) 



if \<k<Nb 
if yt = A^fo 



We recall that £i (Jj.) and £ are function of A and ei {d]^ is the blade efficiency defined in Equations 
2.2 and IJ; its derivative ^^^Jf''^ was already calculated in the Equations 2.17 and 2.19 (Section H]). 



The system of equations p.lq can easily be solved numerically. 



3.2.1 Flat neutron wavelength distribution example 

We take a flat distribution w (A) = jpj^^ between the two wavelengths Ai = lA and A2 = 20A as 
in Section § for the single blade case. In Figure ^the thicknesses of each of the blade coatings and 
each blade efficiency contribution for a 30-layer detector are shown. Three detectors are compared, 
the one of simplest construction is a detector holding 15 identical blades of 0.5 /xm coating thick- 



ness, the second is a detector optimized according to Equation |3.1q for that specific flat distribution 
and the last is a detector that has been optimized for a single neutron wavelength of lOA conform- 



ing to Equations 3.12, 3.13 and 3.14. The fact to have a contribution of wavelengths shorter than 
lOA in the case of the red line makes the coating thicknesses larger compared to the blue curve. 
As a result, frontal layers are slightly more efficient for the optimized detector than for the one op- 
timized for lOA; on the other hand, deep layers lose efficiency. Figure [1^ shows the three detector 
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Figure 9. Thicknesses of the blade coatings (left) and their efficiency contribution (right), for a detector made up of 15 
blades of 0.5 /im, for a detector optimized for the flat distribution of wavelengths and for a detector optimized for lOA. 

efficiencies as a function of neutron wavelength. By comparing red and blue lines, of which the dif- 
ference is plotted on the right plot, the optimized detector gains efficiency on shorter wavelengths 
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Figure 10. Efficiency as a function of neutron wavelength (left) for a detector made up of 15 identical coating thickness 
blades of 0.5 flm, for a detector optimized for the flat distribution of wavelengths and for a detector optimized for lOA. 
Difference between the efficiencies for a detector optimized for a flat distribution and for lOA as a function of neutron 
wavelength (right). 



but loses on longer. Moreover, we notice that the optimization process explained in this section let 
to gain at most 3.5% at short wavelengths while losing less than 1% on longer ones. The weighted 
efficiency over w (A) is shown in Table ^. We can conclude that if we are interested in optimizing a 



Table 3. Averaged efficiency over the flat neutron wavelength distribution (1A-20A) for a detector which contains 
15 identical blades of 0.5 flm, for an optimized multi-layer detector for that specific flat distribution and for a detector 
optimized for 10A(Energy threshold of lOOKeV is applied). 



opt. detect. 


opt. detect, for lOA 


0.5 lJ.m detect. 


0.796 


0.793 


0.764 



detector in a given interval of wavelengths without any preference to any specific neutron energy; 



optimizing according to Equation 3.16 does not give a big improvement in the average efficiency 
compared to optimizing for the neutron wavelength distribution barycenter (about lOA). 
Although the averaged efficiency for the optimized detector in the neutron wavelength range differs 
from the one optimized for lOA only by 0.3% one can be interested to have a better efficiency for 
shorter wavelengths rather than for longer. It is in this case that the optimization process can play 
a significant role. On that purpose let's move to the following example. 

3.2.2 Hyperbolic neutron wavelength distribution example 

We consider a hyperbolic neutron wavelength distribution between Ai = lA and A2 = 20A. 



This optimization aims to give equal importance to bins on a logarithmic wavelength scale. The 

barycenter of the wavelength distribution corresponds to J^^ w{X)XdX = 6.34A. 

In Figure 1 1 are shown the thicknesses of each blade coating and the efficiency as a function of 
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the depth direction in the detector for a 30-layer detector. Five detectors are compared, the one of 
1.2 Hm coating thickness, a detector optimized according to Equation 3.16| for that specific hyper- 
bohc distribution, a detector that has been optimized for a single neutron wavelength of lOA, 1.8 A, 
and for the barycenter of the distribution. 
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Figure 11. Thicknesses of the blades coatings (left) and their efficiency contribution (right), for a detector made up of 
15 identical coating thickness blades of l.ljlm, for a detector optimized for an hyperbolic distribution of wavelengths 
and for a detector optimized for lOA, 6.34A and for for 1.8 A. 
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Figure 12. Efficiency as a function of neutron wavelength (left) for a detector made up of 15 identical coating thickness 
blades of 1 .2 jlm, for a detector optimized for an hyperbolic distribution of wavelengths and for a detector optimized for 
lOA and for 1 .8A. Difference between the efficiencies for a detector optimized for a flat distribution and for lOA, 6.34A 
and for 1.8 A as a function of neutron wavelength (right). 



Table 4. Averaged efficiency over the hyperbolic distribution for a detector which contains 15 identical blades of 
l.lpim, for an optimized multi-layer detector for that specific distribution and for a detector optimized for lOA, 6.34A 
and for 1.8A (Energy threshold of \QQKeV applied). 



opt. detect. 


opt. lOA 


opt. 6. 34 A 


opt. 1.8 A 


l.lpLm detect. 


0.671 


0.641 


0.664 


0.639 


0.597 
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Figure 12 shows the five detector efficiencies as a function of wavelength and their difference on 
the right plot. By comparing the red (optimized detector) and the blue (detector optimized for lOA) 
lines, of which the difference is plotted in blue on the right plot, we notice that the detector opti- 
mized for such a distribution gains about 9% efficiency at short wavelengths and loses about 4% at 
high wavelengths. A detector conceived for short wavelengths, such as the one represented by the 
pink line, has an opposite behavior instead. The distribution optimized detector gains efficiency 
for long wavelengths reaching about 11%. Moreover, a detector optimized for the barycenter of 
the neutron wavelength distribution, instead does not differ more than about 4% over the whole 
wavelength interval, as in the case of a uniform distribution. By only comparing the averaged ef- 
ficiencies, shown in Table 0, it seems that there is not a big improvement in the detector efficiency 
which is only about 3% for both lOA and 1.8A optimized detectors with respect to the distribution 
optimized detector. On the other hand, the optimization procedure, explained in this section, shows 
that it can lead to a significant efficiency improvement in certain neutron wavelength ranges. Fur- 
thermore, as in the case of a flat distribution, a detector optimized for a distribution according to 
Equations 3.16, does not show significant improvement in performances with respect to a detector 
just optimized for its barycenter. 



4. Considerations on solid converter Pulse Height Spectra 

The physical model taken into account in [ [III ] and in [12] can be used as well to derive the analyti 



cal formula for the Pulse Height Spectra (PHS). A similar work was done in [ ]12| ] (see Appendix ^) 
where only Monte Carlo solutions were shown; here we want to use analytic methods to understand 
the structure of the PHS. 

We make approximation mentioned in the introduction and we assume either a simplified stopping 



power function (see Section |4.1.2| ) or one simulated with SRIM [ |18| ] for the neutron capture frag- 
ments. 

We calculate the probabihty for a particle emitted from the conversion point to travel exactly a 
distance L on a straight line towards the escape surface (see Figure |]). This distance L is related 
to the charged particle remaining energy through the primitive function of the stopping power. We 
will demonstrate that even under strong approximations of the stopping power function the model 
still predicts quite well the important physical features of the PHS. 

4.1 Back-scattering mode 

The probability for a neutron to be captured at depth {x,x + dx) in the converter layer and for the 
capture reaction fragment (emitted isotropically in 471^^) to be emitted with an angle (p = arccos(M) 
(between (m, m + du)) is: 

{i • Ze^^'^ dxdu li x<d 
(4.1) 
if x>d 



Where £ is the macroscopic cross-section already defined in Equation and d is the layer thick- 
ness. 
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The fragment will travel a distance L across the converter layer if L = ^ and u G [0, 1]. 

The probability for a particle to travel a distance (L,L + dL) across the layer is then given by: 



P(L)dL= dx du5(--L]p(x,u) = {^^-^^ ) I 



(4.2) 



It is sufficient to replace £ with ^.^^^^ if neutrons hit the layer under the angle 6 with respect to the 
surface (see Figure |l]). The demonstration is equivalent to the one that can found in | U] for the 
efficiency function, in the PHS calculation p{x,u) has to be changed as follows: 

( ^.Ze'^'^^^du ifx<d 
p{x,u,e)dxdu = <^ ™^ ~ (4.3) 

if x>d 

If E{L) is the remaining energy of a particle that has traveled a distance L into the layer, is 
the stopping power or equivalently the Jacobian of the coordinate transformation between L and E. 
Once P{L)dL is known we can calculate Q{E)dE, therefore: 

Q{E)dE = P{L{E)) -j^-dE (4.4) 

I dL I 

where Q{E)dE is the probabihty that an incident neutron will give rise to a release of an energy 
(£■,£■ + dE) in the gas volume; hence it is the analytical expression for the PHS. 

4.1.1 PHS calculation using SRIM output files for Stopping Power 

We take the case of the '"B reaction as example, however results can be applied to any solid neu- 
tron converter. We recall the energies carried for the 94% branching ratio is Eq = lAlOKeV for 
the a-particle and Eq = 830^'^^ for the ^Li; for the 6% branching ratio, Eq = mOKeV for the 
a-particle and Eq = lOlO^eV for the ^Li. Referring to Equation 14, the stopping power used 



here was simulated with SRIM [|T^] (see Figure |^) and L{E) obtained by numerical inversion of 
the stopping power primitive function, i.e. the remaining energy inverse function. 
The complete PHS can be obtained adding the four PHS in the case of '"B according to the branch- 
ing ratio probability: 

Qto,{E)dE = {0.94 ■{QaiE) + QiuiE)) +0.06- {Qa{E) + Q^^iiE))) ■ dE (4.5) 



Consequently the efficiency for a single layer can be calculated by: 

£{ETh)= Qtot{E)dE (4.6) 

J Eti, 



Where Ejh is the energy threshold applied to cut the PHS. This result is fully in agreement with 
what can be calculated by using the Equations in |[TT|]. 

A Multi-Grid-like detector [P was used to collect the data at ILL-CT2 where a monochromatic 
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neutron beam of 2.5A is available. This particular detector has the peculiarity that in each of 
its frames blades of different coating thickness were mounted; as a result the simultaneous PHS 
measurement for different layer thicknesses has been possible. The blades are made up of an 
Aluminium substrate of 0.5mm thickness coated on both sides by an enriched ^*'B4C layer [pj|]. 
Thicknesses available in the detector were: 0.50/xm, 0.75 jXm, 1 jxm, 1.5 jXm, 2}Xm and 2.5 pLm. 
In our calculation we are not taking into account several processes, such as wall effects, gas amplifi- 
cation and fluctuations, space charge effects, etc. but only the neutron conversion and the fragment 
escape. Moreover, while the calculation has an infinite energy precision, this is not the case on a 
direct measurement because many processes give a finite energy resolution. 
In order to be able to compare calculations and measurements, after the PHS were calculated for 
the thicknesses listed above, we convolve them with a gaussian filter of a = \OKeV . The measured 
PHS were normalized to the maximum energy yield {\110KeV). An energy threshold of \%OKeV 
was applied to the calculation to cut the spectrum at low energies at the same level the measured 
PHS was collected. 



We compare calculated and measured PHS in Figure 13; we can conclude that the model gives 
realistic results in sufficient agreement with the experimental ones, to be able to describe its main 
features. 




0.5 um calc. 

0.5 um meas. 

0.75 um calc. 

0.75 um meas. 

1 um calc. 

1 um meas. 

1.5 um calc. 

— 1.5 um meas. 

2 um calc. 

2 um meas. 

2.5 um calc. 

2.5 um meas. 



1000 
E (KeV) 



2000 



Figure 13. Comparison between a PHS calculated and one measured at ILL-CT2 on a 2.5A neutron beam using a 
Multi-Grid like detector [hp where were mounted blades of different thicknesses. 



4.1.2 PHS calculation using a strong approximation 

A fully analytical result that does not appeal to experimental or SRIM-calculated stopping power 
functions, can be useful to understand the PHS structure and to determine its properties. 
The stopping power functions ^ can be approximated by a constant in the case of an a-particle 
and with a linear dependency in L for a ^L/-ion. As a result the energy dependency as a function of 
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the traveled distance L is given by: 



Ea{L) 



-^{L-R) if L<R 
if L>R 



(4.7) 



And equivalently for the ^L/-fragment: 



ELi{L) 



'^{L-RY if L<R 
if L>R 



(4.8) 



Where R is the particle range and £0 its initial energy. 

In Figure [IJ are shown the stopping power functions ^ for '''fi-reaction fragments and their inte- 
gral E{L), in the case of using SRIM (solid lines) and in the case we use the expression displayed in 



the Equations 4.7 and 4^.8 (dashed lines). By substituting Expressions 17 and 4.8 into the Equation 



alpha 1470KeV SRIM 
Li 830KeV SRIM 
alpha 1770KeV SRIM 
Li lOIOKeVSRIM 
alpha 1470KeV approx. 
Li 830KeV approx. 
alpha 1770KeV approx. 
LI 1 01 OKeV approx. 




-alpha 1470KeV SRIM 
-Li 830KeVSRIM 
-alpha 1770KeV SRIM 

Li lOIOKeVSRIM 
-alpha 1470KeV approx. 
-Li 830KeV approx. 
-alpha 1770KeV approx. 
-Li lOIOKeV approx. 




L(um) 



L(um) 



Figure 14. Stopping power ^ and its primitive function E{L) for '"S-reaction fragments, solid curves are the functions 
obtained from SRIM dasiied lines are their approximated behaviors. 



44| we obtain a fully analytical formula for the PHS. It has to be pointed out that each relation, valid 
for L < /?, is valid in the range E < Eq. Hence Equations ^ and [4.10| hold for E < Eq. The two 



formulations in Equation 4.2 for L < J and L> d, translate in two different analytical expressions 
for Q{E) for E < E* and for E > E* , with d = L{E*). 
For the a-particle: 



Q{E)dE= < 



■+d)e'' 



+ R[f 



dE 



E_ 



ifE<Eo{l-i) 
dE if E>Eq{\-^) 



(4.9) 



Where the relation E* 



is derived from d = L{E*). 
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For the Li: 



Q{E)dE- 



[l-{{ + d)e-'^'^)dE if £< £0(1- I)' 



(4.10) 



-ZR{ 



*j if E>Eo{\-{f 



Where, again, the relation E* = Eq{\ 



is derived from the condition d = L{E*). 



Figures |15J, |16| and |17| show the calculated PHS obtained by using the SRIM stopping power func- 



tions and the approximated one displayed in the Expression [1^.7| and [4.8|for 0.2 pim, 1 pim and 4 pim 
respectively, when neutrons hit at 90° the surface and their wavelength is 1.8 A. They show similar 
shapes that differ in some points; e.g. focusing on the 1410 KeV a-particle, the fact that the approx- 
imated E{L) function (see Figure [I^ differs from the SRIM one at high L leads to a disappearance 



of the PHS rise at low energies; it is clearly visible in Figure [17. 

We see that as d increases what looked like a single peak splits into two peaks. While one peak 
stays constant at the highest fragment energy Eq the second one moves toward lower energies when 
the layer thickness increases. This is important when trying to improve the neutron to gamma-rays 
discrimination by creating a valley that separates them in amplitude. 
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Figure 15. Calculated PHS using SRIM (left) and approximated (right) stopping power functions for a single back- 
scattering layer of 0.2/im for 1.8 A and 90° incidence. 



In order to understand the PHS structure, we define the PHS variable space: on the abscissa axis is 
plotted u = cos((p), where (p is the angle the fragment has been emitted with, and, on the ordinates 
axis, is plotted the neutron absorption depth x. u £ [0,1]; x € [0,d] if d < R or x £ [0,R] if d > R 
because a neutron can only be converted inside the layer and, on the other hand, if a neutron is 
converted too deep into the layer, i.e. x > Rno fragments can escape whatever the emission angle 



would be. In Figure |18|, on its left, the variable space is shown; an event in the A-position would 
be a fragment that was generated by a neutron converted at the surface of the layer and escapes 
the layer at grazing angle. An event in the position B represents a fragment that escapes orthog- 
onally the surface and its neutron was converted at the surface. An event in C means a neutron 
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Figure 16. Calculated PHS using SRIM (left) and approximated (right) stopping power functions for a single back- 
scattering layer of 1 flm for 1.8 A and 90° incidence. 
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Figure 17. Calculated PHS using SRIM (left) and approximated (right) stopping power functions for a single back- 
scattering layer of 4 jxm for 1.8 A and 90° incidence. 



converted deep into the layer with an orthogonal escaping fragment. The straight line x = L{E)-u 
characterizes the events with identical escape energy E, that contribute to the same bin in the PHS. 
The straight line characterized by a; = /? • m is the horizon for the particles that can escape the layer 
and release some energy in the gas volume. To be more precise events that give rise to the zero 
energy part of the PHS lie exactly on the line ^ = L{E = 0) =R because they have traveled exactly 
a distance R in the converter material. On the other hand, events that yield almost the full particle 
energy £"0, will lie on the line identified by = L{E = Eq) = 0. 

The events that generate the PHS have access to a region, on the variable space, identified by a 
triangle below the straight line x = R-u (see Figure Q. 



lid > R the variable x can explore the interval x G [0,R]. This is the case of the PHS in Figure |17| , 
where d = Ajim, Ru{S3QKeV) = 1-7 /im and Ra{uiQKeV) = ^-^IJ-m. We take the two particles of the 
94% branching ratio of reaction as example. 



If d <R, the variable x can explore the interval x G [0,d] (see Figure 19), thus the domain is now 
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Figure 18. PHS variable space sketch (left) and PHS variable space in the case d >R (right). 




a trapezoid. The events near the line ^ = L(E*) = d, which is the switching condition found in 
the Equations and 4. lO] , give rise to a peak because this line has the maximum length available. 
Thus, we expect a peak in the PHS around E* . This is shown in Figure 16 where d = \ jim and, 
again, Ru(?.mKeV) = 1-7 pim and Ra(\MGKey) = 3.4/im, the peaks that originate from the condition 
^ = L{E*) = d for the two particles are indicated by the arrows. 

\f d >R, the peak occurs for L{E*) = R, that is, zero energy. This is problematic for /-ray to 
neutron discrimination. 

If d << R, the variable space is compressed and the straight lines identified by ^ = L{E*) = d 
and ^ = L{E) = become more and more similar. The two peaks approach and, the more d is 
negligible compared to R, the more the two peaks appear as one single peak (see Figure 15). 
If we want to avoid a strong presence of neutrons in the low energy range of the PHS, where we 
know the y-rays contamination is strong, it is important to try to get the second peak higher than 
the energy threshold {Eth)- This implies that the thickness d of any layer in the detector should 
obey d < L{Ejh) for the L corresponding to the particle with the smallest range. This can be a con- 
tradictory requirement with efficiency optimization in which case a compromise between y-rays 
rejection and efficiency has to be found. 



4.2 Transmission mode 

Equations for transmission mode can be calculated in the same way they have been determined 
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for back-scattering mode by substituting x with d — y m the Expression and x with y in the 
expression 5 — L) . As a resuh we obtain: 



P{L)dL = l^^^ ^ ^' > - (4.11) 

\^ {y-^'^ + {d - {)) dL if L>d 

Hence, Q{E)dE can be calculated as shown already in Section p~T| . 

However the same conclusions can be drawn concerning the qualitative aspects of the PHS, espe- 
cially the position of the two peaks. 



5. Conclusions 



We demonstrated that the sputtering technique, suited to make blades with equal coating thick- 
ness on both sides of the substrate, is well-adapted to make efficiency optimized blades when the 
substrate effect can be neglected for any given neutron wavelength and for any incidence angle 
distribution. Moreover, this result is also valid for a multi-layer detector where several blades are 
arranged in cascade. 

Analytical formulae have been derived in order to optimize the coating thicknesses of blades in 
single-blade and in multi-layer detectors. 

The blade-by-blade optimization in the case of a multi-layer detector for a single neutron wave- 
length can achieve a few percent more efficiency over the same blade optimization but this can lead 
to several blades less in the detector. Moreover, in the case of a distribution of wavelengths, the 
suited optimization from a distribution does not give important improvements in the overall effi- 
ciency compared with a monochromatic optimization done for the barycenter of the distribution. 
On the other hand, the optimization of the efficiency for a neutron wavelength distribution is often 
more balanced between short and long wavelengths than the barycenter optimization. 
We have demonstrated that for our model the analytical expression for the PHS is in a good agree- 
ment with measurements. Moreover, thanks to this model, we understood the overall shape of the 
PHS which can be important if one wants to improve the /-ray to neutron discrimination in neutron 
detectors. 



6. Outlook 

Even though the substrate effect can be neglected in most cases when dealing with a small amount 
of blades, its effect in a multi-layer detector can strongly differ from the results obtained for the 
ideal case of completely transparent substrate. A further step is to take its effect into account. 

A. Formulae in []TI|] 



The relations between the formulae in [ ]11[ ] and the expression used in this paper are the following: 

• the particles ranges L are denoted by R; 

• the branching ratios of the ^^B reaction (expressed by Fp) are Fi = 0.94 and F2 = 0.06; 
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the thickness of the layer is d = Dp. 



Hence, the relation between the expressions y..\\ |2.2|, 2.5 and the formulae in 1 1 1] is 



EridT) = 0.94 • eT{R\^'^%Rf"'")+0.06 ■ SriRf^Rf') = 

(A.l) 

= Si{Df,Li, 0.94) + Si (Df , L2 , 0.94) + S2 (Df , Li , 0.06) + 52 (Df , L2 , 0.06) 
Valid for both equations (18a) (Dp < Li) and (18&) (Dp > L,) in Section 4.2 of [llT]]. In the case 



of the back-scattering mode, equations (25a) and {25b) in []11[], we consider one layer of converter 
and we replace Eridr) into SBsi^Bs) in the expression for back-scattering. 
In a different way, for both equations (18a) and (18^?), we can also write: 

Fp-eT=Sp{DF,Li,Fp)+SpiDF,L2,Fp) (A.2) 

B. Formulae in 



The relations between the formulae in [ ]12p and the expression used in this paper are the following: 

• the macroscopic cross-section (E) is expressed in terms of mean free path / = ^ ; 

• the variable u is denoted by its cosine u = cos(0); 

The formulae (4) in [|l^] corresponds to the Equation p~T| , unless a factor ^, where / = ^■ 



The formulae (3) in [ jl^ ] corresponds to the Equation ^ 
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